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1 Generalities 

Inflationary magnetogenesis has long been recognised as one of the most promising frame- 
works for the generation of today's observed cosmological magnetic fields, see [1-5] , and [6- 
12]. Despite the initial hindrance represented by the conformal conspiracy of conformally in- 
variant Maxwell electromagnetism (EM) and conformally flat Friedmann-Lemaitre- Robertson- 
Walker (FLRW) universes, abundant variations on this theme have been concocted by theo- 
rists, some of the simplest options revealing themselves as among the most elegantly successful 
in this competition. 

So, why inflation? In our opinion, besides the obvious ease with which such models 
can be studied - in their most humble guise all one needs to do is to solve a homogeneous 
second order differential equation - the reason is that the property which is perhaps the most 
intricate to frame within a consistent theoretical model is their pervasively wide coherence 
length, which can stretch well beyond the Mpc mark: cosmological inflation, unlike later 
time causal mechanisms, is automatically equipped with this feature. 

Inflationary magnetogenesis is not without its list of drawbacks: many simple options 
suffer from a backreaction issue, see [13-16]. Very briefly, as inflation amplifies EM quantum 
fluctuations, the overall energy density of the latter ends up being comparable with that of 
the former, thereby implying the breakdown of perturbation theory. When this occurs either 
inflation stops - and if this happens too early then it does not generate enough entropy 
and curvature perturbations to explain today's Universe, or the magnetic fields cease to be 
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amplified and remain too small today 1 ; the possibly simplest viable option [18] appears to be 
successful only in the strong coupling (hence non perturbatively trustable) regime [13, 14]; 
furthermore, the counterproposal of [19] also would not step down this cross [20]. 

In this battle for inflationary magnetogenesis there is, however, a survivor: this is 
the three- form driven inflation [21-23], which can be coupled in an elegant way to EM to 
amplify its initially tiny quantum fluctuations to sizeable amplitudes today, and at scales 
which easily exceed the Mpc scale [24]. In one line, the reason this particular model sails 
the winds where others float or sink, is that within it it is possible to carefully select which 
band of EM modes are to be amplified (for instance, the low momentum ones), thereby 
compressing all available energy density before backreaction into a handful of scales only. 
This can be understood since the vector degrees of freedom in the three-form are excited by 
the coupling (consequently breaking down the duality with a scalar field inflation), resulting 
in new modes with a nontrivial dispersion relation. In particular, there appears a new mode 
that is exponentially growing over a certain range of wavemodes. 

This is, however, not the end of the story. Going beyond background dynamics, where 
only the EM field is treated as a small perturbation, a first order in the metric perturbations 
analysis was developed in [25]. Moreover, such coupling between the background scalar and 
EM can genuinely lead to non-Gaussianities in the Cosmic Microwave Background spectrum, 
see [20, 26-28] and correlations between the scalar perturbations and the magnetic field [19, 
29] . The EM field generated during inflation contributes in general an anisotropic term in the 
energy-momentum tensor, which marks the difference between the two scalar gravitational 
potentials (in longitudinal gauge), and has to be taken into account. The conclusion is that 
such piece partly transmits to a constant curvature perturbation in the following radiation 
era: if this new term, whose absolute value depends on today's magnetic fields, becomes 
comparable to the background, it signals the breakdown of perturbation theory. 

The scope of this work is to revisit this crucial issue for the particular case of the three 
form background. We first review "three-magnetogenesis" in Sec. 2, and the appearance 
of the constant mode in Sec. 3; we apply this analysis to three-form inflation in 4; with 
the tool-set at hand we briefly discuss non-Gaussianities generated in this scheme (Sec. 5); 
conclusions make up the closing Sec. 6. 

2 Three-magnetogenesis 

We review and expand on the model first introduced in [24]; in particular, we complete that 
work by presenting the derivation of the equations of motion in some detail. We also discuss 
alternative prescriptions for the initial conditions of the system. 

2.1 Dynamics 

The dynamics and stability of three-form cosmology have been studied in [23, 30-32] and 
considered also in the cases of nonminimal couplings to gravity [21, 33, 34] or to matter [35]. 
In the following we will focus only on what is new here, that is the coupling to the electromag- 
netic field and the ensuing dynamics of the vector perturbations. Trivially, the three- form 
field can be equivalently described by the dual vector field, and this sense the models are 
examples of vector inflation [36-39]. The equivalent vector description would have a non- 
canonical form, but such vector models have been also studied recently [40-42]. 

1 Backreactions were first mentioned in an earlier work by Campanelli [17], where they were found to be 
unimportant thanks to a clever, albeit entirely arbitrary UV cutoff. 
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Let us first specify the Lagrangian we study. Let be the photon vector potential and 
B tJ,up the three-form. The canonical Lagrangian including both fields is 2 

Ca + C b = -\f\A) - - V{B 2 ) , (2.1) 

where the Faraday forms are computed from an n-form potential N as F(N)p 1 ^^ n+1 = 
(n + l)!<9[ m iV M2 ^ n+1 ]. The components of the dual of the three-form are [43] 

B a = l -e a p l5 B^ s . (2.2) 

The most general Lorentz-invariant, quadratic, second order and £7(1) invariant coupling of 
the two fields is of the form 

Cab = - l -\iF^ v {A)F^{B) (2.3) 
+\ 2 {V^){V U B U ) + (A 3 + \aR)A^ + X 5 A^R^B U + \ 6 A»R^ pa B^ a . 

Here Rp Vpa is the Riemann tensor, and R^ v = R 9 p,p V and R = R^ its Ricci contractions. 
Only the first term respects the £7(1) symmetry of the photon. The second term breaks it 
only partially, leaving the reduced gauge invariance under the transformation A — > A + d<p, 
where <f> is a harmonic scalar field obeying n<fi = 0. If A3 = 0, the remaining three terms 
vanish in flat space, and thus the standard predictions of QED are recovered regardless of 
the dynamics of the three-form. 

We choose to work with only the £7(1) gauge invariant term, and set all Aj to zero, 
except for Ai. We expand the spatial part of the vector potential A^ = (Ao,Ai) in terms 
of its transverse and longitudinal components as Ai = Af + d{A L , where V • A T = 0. The 
three- form is similarly decomposed Bi = Bj + diB L . Since the physical photons correspond 
to the transverse degrees of freedom, we are particularly interested in the vector perturba- 
tions. The line element of the Friedmann universe including rotational perturbations can be 
parametrised as 

ds 2 = a 2 (77)[dr7 2 + c T • dxdr? - (dx • Vd T ) • dx] . (2.4) 

We will work in terms of C T = — c T + d T> , which includes the two gauge- invariant rota- 
tional degree of freedom of the metric perturbations. The vector part of the Einstein-Hilbert 
Lagrangian coupled to the three-form is 



_i M 2 k 2 qT2 v x ^ bT YC ?y 

2 X 



(2.5) 



where we have used the background Friedmann equation and X is the background value of 
the three-form field X = V B 2 in the notation of Ref. [30]. 

Let us now redefine Ai = A. The relevant part of the action (2.5) simplifies to 

S A + S AB = ^J { A t (-d 2 + A) (A T + 2AB T ) - A A [a - 2A L ' + 2A (b - B L '^j 

-A L 'A [A L ' - 2A (b' - £ L ')] } a 4 (??)d 4 x . (2.6) 



2 In terms of the dual vector, V = V(6B 2 ) and F 2 {B) ~ (V • B) 2 
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Varying with respect to Aq, we get Ao = A L ' — X(Bq — B L '), and the action further reduces 
to 

C A+AB = \ [A T (-d 2 + A) (A T + 2AB T ) + A 2 (b - B L ') A [b, - B L ')] , (2.7) 

where rj is conformal time, and A = 5 lJ didj. The scalar polarisations of the photon are 
nondynamical, as expected. Notice however that while the scalar part of the interaction 
decouples from the EM field as it should, it nontrivially contributes to the effective speed of 
sound of the three-form. 

The equations of motion for the vector degrees of freedom become (omitting from now 
on the superscripts T, since all vectors are considered to be transverse) 

M^k 2 C + 2V x (B- XC) = 0, (2.8) 
Vx(B - XC) - \X(-d 2 + A)A = , (2.9) 
(-<9 2 + A)(A + 2AB) = 0. (2.10) 

We note that in the absence of the coupling, both the metric and three-form vector perturba- 
tions are nondynamical. We can eliminate C and obtain a closed equation for the three-form 
perturbation in Fourier space: 



B = 0, (2.11) 



/ = »*£(^g-i|. (2.12) 



where 

f — 9 \ 2 

V x V M*k 2 

and B is the Fourier transform of B and similarly we'll define A, 

A(K V ) = j^j- 2 J d 3 xA(x, n)e-^ , B(k, V ) = ^_ J d 3 XjB(X; ^-a** . (2 . 13) 

Thus the three-form rotational modes propagate with a nontrivial dispersion relation; they 
can therefore in principle be significant even at large scales. 

We are in the position to write an equation for the EM potential and the three-form 
only as A" + k 2 A = F(B), where F(B) = -2A(<9 2 + k 2 )B: the EM potential behaves as 
a harmonic oscillator driven by an external force. Moreover, we can finally present the 
autonomous evolution equation for A, 

(-d 2 + A) [1 - f(-d 2 + A)] A = . (2.14) 

Each transverse degree of freedom in momentum space obeys 

A^ + 2^A^ + y(f" ~ 1 + 2fc 2 /K' + 2^14' + ^(^^ + k 2 )A = . (2.15) 

Fourth order time derivative appear in our equations, which may signal the presence 
of ghosts and/or negative norm states. Since this is an effective field theory, and as such is 
only valid at low energies k <C V^ 4 (Vq being the UV mass scale of the model), as long as 
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we confine ourselves to such energy range we are allowed to not worry about this potential 
drawback. 

In order to be able to follow the analysis analytically, we consider de Sitter solutions with 
constant comoving field X. The two classes of fixed points [30] are named B (corresponding 
to X 2 = ±|Mp), and C (corresponding to the minima of the potential). Depending on the 
shape of the potential these can be stable or unstable; for this discussion what we need to 
know is that in both cases Eq. (2.15) much reduces to the simpler 

A^ -{j o - 2^ 2 ) A" ~ k 2 (j - fc 2 ) A = , (2.16) 

/o is given by Eq. (2.12) evaluated at the fixed point. We can readily solve this equation to 
obtain 

A(r)) = A! cos (krj) + A 2 sin {kr,) + A 3 e rkri + A A e~ Tkr! , (2.17) 

where Aj are constant vectors and T 2 + l = l/(/o/c 2 ). Hence, we have two additional solutions 
(as the equation is fourth order) which can be oscillatory if T 2 < 0, or, more interestingly, 
exponentially growing / decaying when T 2 > 0. 

2.2 Examples 

A perfect example of the features of this model is given by the exponential potential V = 
Vbexp (— £X 2 /Mp). At the fixed point B, which is stable as long as £ is positive, we have 

K 2 - K 2 

r 2 = . K 2 \ 2 2 , (2-18) 

where we have defined A 2 = 8A 2 /3/(l - 8A 2 /3) ~ 8A 2 /3 and k\ = 8£V/(3M 2 A 2 ) ~ 
£V/(\ 2 Mp). Only for k < k\ there exist an exponentially growing solution: k\ acts in 
this case at the de facto UV cutoff of the theory. We have finally defined k = k/T-L e for con- 
venience, where 7i e is the highest mode accessible by inflation - we further denote % = a' /a, 
where a is the scale factor of the Universe and the prime stands for conformal time derivative. 
The second fixed point C can be an attractor if £ is negative. In this case one sees that 

rt-^-^. (2.0) 

Thus, when this point is an attractor, there is no instability of the vector modes. If we are 
inflating at a saddle point (i.e., £ > 0) vector modes will be unstable at scales k 2 < A f . 

The properties of the fixed points for many classes of potentials and their T 2 + 1 are 
summarised in tables 1 and 2. 

One can therefore quite straightforwardly realise when the instability for the vector 
potential A, which is parametrised by T 2 , is present, and where. If large scales, k — > 0, are 
exponentially amplified, then AX 2 X 2 /Mp < 1 or else we only have plane waves. Also at large 
momenta we recover plane waves because T 2 — > —1. Hence, there is in general a low energy 
band which is unstable given the above condition. 

In addition to the turning point where T 2 changes sign, one should worry about one 
more singularity in T 2 which appears for the mode k which nullifies its denominator. The 
singularity could lie within the exponentially growing band, which is problematic since it 
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Table 1. Stability of the de Sitter fixed points of type B and the characteristic width of instability 
of the vector potential in four classes of models. U - unstable; S - stable. 
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Table 2. Stability of the de Sitter fixed points of type C and the characteristic width of instability 
of the vector potential in four classes of models. U - unstable; S - stable. 



would mean that to that mode corresponds infinite energy density. One can preventively 
understand whether this will be the case by looking at the sign of V t xX: if this is positive 
the divergence will appear; in the opposite case VxX < the troublesome scale is beyond 
the UV cutoff. This uncontrolled instability is consistently steered clear of in what follows. 

We choose then to work with the exponential potential V = Vq exp(— £X 2 /Mp) because 
it does not require a very high degree of fine-tuning; a similar phenomenology arises in the 
Ginzburg-Landau case V(X) = V (X 2 - C 2 ) 2 /Mp, but only at the fixed point B. 

2.3 Some results 

The action (2.7) is canonical for the EM field; thus, since we regard the coupling to the 
three-form as a perturbation, we are allowed to proceed by quantising canonically [44]. The 
solutions of interest are those which grow with conformal time, so we will choose A\ = Ai = 
A4 = - this might appear contrived, but it is clear that, were these solutions kept, they 
would very rapidly loose any relevance since, in comparison to the growing ^3 piece, A\ and 
A2 are just plane waves, while Aa decays exponentially. 

At this point we need to fix the initial conditions (IC) of the dynamics for the EM vector 
potential A; we are presented with two possibilities: either we set them at some initial time 
(e.g. the beginning of inflation r/j), for all modes, or we normalise each mode to the Bunch- 
Davies vacuum y/2kA = e~ lkri at horizon crossing. This last choice corresponds to a coupling 
which dynamically becomes relevant at some time, while being inefficient for smaller scales. 
This option makes it clear that for all modes the specific IC do not matter, while they can 
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be consistently chosen to be in vacuum at any early time without affecting the final result. 
However, this demands a more elaborate theoretical motivation for the time-dependence in 
A. 

On the other hand, the former possibility is simpler, but perhaps not the most physical 
nor appealing one, for it introduces an explicit and very strong dependence on the IC of 
inflation, which is antithetic to what inflation is brought in for. Moreover, there remains the 
formal question on how to consistently choose a vacuum state since the effective coupling is 
never turned off. This, actually, can be fixed if we work in the small coupling regime: in this 
case it is consistent to choose the Bunch-Davies solution as the approximate IC. 

One final consideration concerning the IC for the system. We work with the simplifying 
assumption that X is a constant, which implies that / is so. Then, the equation of motion 
for B (2.11) does not admit plane wave solutions at any time. Nonetheless, it is quite easy 
to convince ourselves that in more realistic models, when fk 2 — > ±00 as kr\ — > —00 (which 
depends on the shape of the potential) then the rotational perturbation B can, and in fact 
must, be initially in its Bunch-Davies vacuum, which in turn is realised for A. That is, the 
exponent V in reality is time-dependent as well, and flips from imaginary to real at a given 
time; we take the latter to be rji for simplicity. 

IC at horizon exit - Let us work out the last case first. The IC fix the value of ^3 
and give us the following solution to work with: 

A = -^=e v ^ +l) IC at r, = -l/k .; (2.20) 



Before moving on to the actual discussion, the key definitions we will need are the magnetic 
power spectrum 



where a is the scale factor of the Universe, and the energy density in electromagnetic field: 

^ = ^/f* 3 ^ + *V| 2 ], (2-22) 

where the first term is related to the electric field, and the second one is the magnetic energy 
density proper. Recall that for a de Sitter solution conformal time spans the range [— 00, 0], 
and % = — l/r?. We can immediately write down the corresponding expressions for the energy 
density associated with the EM field at the end of inflation as 

Ut J Akk 3 (r 2 + 1) e 2 ^ 1 -*) . (2.23) 



P cm - - 9 a 4 



8vr 2 



where of course K e = 1 by definition, and where we can clearly ignore the second piece in 
the exponent, being greatly smaller than the first, positive, contribution. Also, we work with 
small A, which means A <C 1 and r> 1; this implies that in general E 2 3> B 2 . 

The energy density Eq. (2.23) is expressible in closed form. The limits of integration 
are k = Ki = exp(— N), where N is the total number of e- folds inflation lasts for, and 
k = ka, beyond which the solution (2.20) does not apply any more. The final expression to 
be evaluated at those limits, reads 
-1/4 4 ( 2r 

Pem = 8^4 f { TTT 2 (1 + A ' + A ') + ic * t 1 + A ' + (! " *) A 1 ^ ( 2r " 2 *) 

[l + A 2 + (l + i)A 4 ] Ei(2T + 2i)} . (2.24) 



-te 



-2i 
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This expression is not particularly illuminating in this form, but can be expanded for small 
A to yield 



Pem = -^ K i A 3 e 2/A = ^L A 3 e 2/A (2>25) 

Notice how the cutoff scale k\ appears, where normally one has ~H e , and how the exponent is 
inversely proportional to A. In the same fashion, we can express the spectrum as seen today, 
diluted by a factor (a e /ao) 4 , as 

,2,0 _ k*_ 2T _ 2/VA 2 +(fe/fe A ) 2 (o 26) 

In order to have a low-scale cutoff the scale of the three-form potential needs to be low. 
Moreover, since T oc 1/A we want a small coupling constant for the most efficient amplifica- 
tion. 

In Fig. 1, left panel, we present the regions in the parameter space (A, &a) for which 
the ratio p em /px stays below one at the end of inflation (blue region) and goes beyond one, 
thereby leading to instability in the de Sitter background (red region). Superimposed are the 
contours for a given value of the magnetic power at 1 /Mpc today, where the regions for which 
$B — 10~ 15 Gauss and 5g > 10~ 25 Gauss are in evidence (green and yellow, respectively). 
These are values that correspond to the magnetic fields observed in the intergalactic medium, 
and a typical value for a successful seed to be fed the magnetohydrodynamical plasma at late 
times. The right panel of Fig. 1 plots the power spectrum today; notice again the knee at high 
energy and the following very rapid decay. Notice how such a spectrum also nicely helps in 
avoiding the stringent constraints which late time cosmology poses on smaller scale fields [45] . 
These figures were obtained for around N = 66 e-foldings of inflation, % e /a e ~ 10 13 GeV, 
and assuming the Universe was always dominated by radiation from the end of inflation 
onwards, but they are by construction fairly independent on these specifics. 

IC at the beginning of inflation - For clarity, we discuss separately the calculation 
with the alternative prescription for the initial conditions. Let us set: 

A = -^e rk(r i-^ IC at n = r)i . (2.27) 
V2k 



The energy density is then computed from the integral 

U 



J dKK 3 (r 2 + 1) e 2r«(V«i-i) . (2.28) 



The integral (2.28) instead is not simply dealt with, and we must resort to numeric methods 
to solve it. In Fig. 2, left, we present the same regions in the parameter space (A, k\) for 
exp(— iX 2 ). Fig. 2, right, plots the power spectrum today, for A = 10~ 2 and k\/k m \ n = 10 2 . 

As one can see by looking at these figures, and as it can be expected based on Eq. (2.27), 
there is some fine tuning in this case; there is a very strong sensitivity (doubly exponentially 
in terms of the total number of e-folds N) to the initial time, and only scales not too far 
from the first ones available to inflation can be amplified efficiently without leading to a 
catastrophic backreaction effect. We picked the representative scale k = 3/Gpc, and the 
largest mode is only around 1/Ho ~ 4Gpc. By looking at the spectrum we realise how only 
a handful (not more than a couple, in this case) of orders of magnitude in k will be amplified 
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log e B e 1 / MpC 



2 

5 




Figure 1. Le/£. Allowed (in light blue) and disallowed (in light red) regions in the (log A, log (k\Mpc)) 
parameters space; in the green region the generated magnetic field can explain magnetisation at 
and beyond the Mpc scale, while in the yellow region some degree of dynamo would be required. 
Superimposed are contours of constant log 8^ today, at 1/Mpc (in Gauss units). Right. The power 
spectrum today, 5^ (Gauss), for A = 10~ 2 and — 10 4 /Mpc. 




Figure 2. Left. Allowed (in light blue) and disallowed (in light red) regions in the 
(log A, log (fcA/^min)) parameters space; in the green region the generated magnetic field can explain 
magnetisation at and beyond the Mpc scale, while in the yellow region some degree of dynamo would 
be required. Superimposed are contours of constant logS% today, at 3/Gpc (in Gauss units). Right. 
The power spectrum today, Sg (Gauss), for A = 10 -2 and fcA/fc m ; n = 10 2 . 



efficiently; for instance, the strength at 10/Gpc is 10~ 21 Gauss, but drops to 10 Gauss at 
just 30/Gpc. The scales we amplify therefore heavily rely on what the largest scale is. With 
the alternative prescription of a time-dependent coupling the sensitivity to initial conditions 



- 9 - 



can in principle be eliminated. Notice that here inflation lasts about 60 Hubble times. 

From the allowed regions in Fig. 2 it seems that we could push the coupling A ~ A 
to arbitrarily small values, and still be able to obtain sizeable power at the scales we need. 
This is only true in part: since A also enters the UV cutoff, and as A — > one needs to 
keep A 2 A;^ — £Voe _2 £/ 3 /Mp fixed, the window of modes which can be boosted has to be 
kept narrow. Lastly, one could consider models with either very high or very low potential 
slope parameters £: this would help reconciling suitable UV cutoffs with weaker interaction 
strengths A. 

Having now shown that the three-form magnetogenesis can avoid the backreaction prob- 
lem at the zeroth order, we turn to consider the linear perturbations more carefully. 

3 The framework for perturbations 

We now consider the scalar sector of metric perturbations, and review the results of [25] 
following closely their presentation. 

3.1 Metric perturbations and anisotropies 

The starting point of the analysis of [25] is the perturbed metric in longitudinal gauge: 

ds 2 = a 2 [(1 + 2$)d?? 2 - (1 - 2f)dx 2 ] . (3.1) 

In order to obtain the perturbed Einstein field equations SG^ = 8irG5T^ u we need an expres- 
sion for the energy- momentum tensor (the Planck mass is 8ttG = Mp); we can decompose it 
with respect to the velocity of the FLRW background, = a _1 (l — 3>, 0), and to the metric 
g^ u as is customary as: 

TZ = (Pern + PemK^ " PemiT + «"C + + > (3-2) 

where we take the EM field to be half order in perturbation theory, where, thence, only the 
unperturbed parts of the velocity and metric g^ u will appear. 

In standard inflation the background fluid is described by a scalar field, p, expanded 
as (p(ri, x) = <p(r)) + Stp(r],x.), with a standard kinetic term and potential. In this case, the 
background equations of our system are 

4vrG^' 2 = U 2 -W , (3.3) 
0" + 2U<p' + a 2 V tV = . (3.4) 

It is straightforward to write down the perturbed Einstein equations now, which, in 
Fourier space and in the presence of an electromagnetic field are: 

3W + {2U 2 + H')<S> + = -4irG(ip'6<p' + V tV a 2 S(p) - 4nGa 2 Pem , (3.5) 



tf" + 2W + H& + {2U 2 + H')& - y ($ - *) = 4ttG(VV - V,ya 2 6tp) + 4irGa 2 Pcui , 

(3.6) 

+ %<$> = 47tG<p'5ip - AirGa i p q cmj , (3.7) 
k 2 {<5>- #) = -87rGa 2 n 5 . (3.8) 
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We define p em (\t), p e m(k) and q era j(k) as the EM field energy density, pressure and Poynting 
vector (all these in momentum space). We also define Ils(k) = —3kkjHj mi (k)/2 as in [25]. 
Notice that the background equation (3.3) has been used here. The anisotropic contribution 
Us in Eq. (3.8) is what prevents the usual identification $ = ^. 

Eliminating V tip through the background equation (3.4), using Eq. (3.7) to obtain <p'6(p, 
and finally adding (3.5) and (3.6) we derive an equation (Bardeen) for ^ alone as 

2H<p" 



V + 2 [H - 



1_ 
0' 



0' 



(3.9) 



where every <I> instance has been replaced by VP through Eq. (3. 
collects all contributions from the EM field, and reads 



a^k* 



- 7,(pcm ~ Pcm) + [ 2U + 



2 W 



i Qem j 



3 



Cp' ) k 2 a 



The source term 5 em 
k 2 

(3.10) 



Since we are interested in magnetic fields which today have extremely large wavelengths, 
we focus our attention to scales for which x = k\rj\ <C 1 (the conformal time in almost de 
Sitter background being negative) and we develop everything in the slow roll approximation; 
the slow roll parameters themselves are defined as 

n 2 -n' = eiH 2 and e[ = 2e 1 (2e 1 + 3e 2 )n . (3.11) 

These can be recast in another useful form for our subsequent computations as 

n' = U 2 {1 - ei) and U" = 2U Z {\ - 2e x - ej - 3eie 2 ) . (3.12) 

Time, scale factor, and the Hubble parameter are thence related as a oc l/r/ 1+ei and % ~ 
-(1 + ei)h- 

Putting all this together, we arrive at the master equation 
j2, Tl o/, , q, ^,t, r 2(2ei + 3e 2 )" 



d^* 2(ei + 3e 2 ) dtf 
+ — L — + 



dx 2 



dx 



1 



x* 



where 



S ~ 



87rGa 2 ^ 2 



n</ k 2 k 2 

_A_ + 2 (1 - 2ei - 3e 2 ) U s - —U s - —{p cm - P em ) + 3fe 2 Q e 



where we defined 



Qem — 



iWeig, 
k 2 a 



em j 



(3.13) 



(3.14) 



(3.15) 



and iS em = S era /k 2 . We have expanded everything to order x 2 inside the parentheses, but 
have kept the the first nonzero slow roll order only for the term 11$ /x 2 , the reason for which 
will become clear below - this is a consistent expansion for our purposes. 

It is possible to show [25] that the flow Q is related to the other quantities in the source 

by 



Q'er 



,-n s - p« 



H 3" 

Later we will also use the fact that, of course, P e] 



m 2 - W 

U 2 



Qe 



(3.16) 



/3. 
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3.2 The source unveiled 

In order to solve Eq. (3.14) one needs to compute Us; this can be done if the time evolution 
of the source generating the anisotropy is known, which in turn, in the EM case, can be 
read off the solutions of the equation of motion for the EM vector potential coupled to the 
inflaton. The simplest, and analytically solvable, model in this sense is given by 

C = -\f(v)F^F, u , (3.17) 

In Coulomb gauge Aq = 0, djA^ = it is straightforward to solve the Euler-Lagrange 
equations for the Fourier modes 3 fA 

(fA)"+(k 2 -^jfA = 0, (3.18) 

where the coupling function f((p) is parametrised as a power law as: 

f(v)ocrf. (3.19) 

Typically the value of 7 is restricted to be within [—2, 2], both because most quantities 
(power spectra, energy densities, etc.) are well-defined in this range, and the EM field 
remains subdominant during sufficiently long-lasting inflation. In fact, as for the latter, it 
can be shown that slightly larger values of I7I < 2.2 are allowed (see for instance [13]). 

Skipping the details of quantisation, the solutions to (3.18) are given by 

A(k, n) = \J\ [Ci(7) Vi/2(z) + C 2 (7)J- 7+ i/2(aO] , (3.20) 
where J v are Bessel function of order u, 

f^ e -mi/2 e i7r( 7 +l)/2 

d( 7 ) = J-r , and C 2 ( 7 ) = J-r • (3.21) 

V 4 COS7T 7 V 4 COS7T 7 

The long wave expansion for this solution reads 

A(k,rj) ~ ^[c 1 (7)^ + di(7)x 7+2 + c 2 (7)x 1 - 7 + d 2 (7)^ 3 ~i , (3.22) 
V k L J 



with 



c l(7 )= e 7 /2 , ^ , 4(7) = -^, (3.23) 

1KV cos(^7) 2 7-ir(7 + l/2) 7 + 1/2' 1 ; 

i7r( 7 +l)/2 /^M c („) 

C2 ( 7 ) = v^/ 4 d{) = _^vtj_ (3 _ 24) 

2W; cos(^ 7 ) 22-TT(3/2- 7 ) 3/2 - 7 1 ' 
Every term of the source S em can be extracted from the general expression 

V^/'fv*-^) , (3-25) 



3 Notice that in order to apply our formulae in the previous and following sections we need to canonically 
normalise the field. 
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At this point, overflying the details of the hard part of the calculation, three steps are 
necessary in order to derive a useful expression for the source term. First, since it can be 
shown that only the time dependence of this term is relevant for us, the most direct way to 
do so is to compute the power spectrum 



<0|n s T( q ,, 7 )n s (k ) r/)|0> = (27r) J P n (fe,r?)5(q-k) 



(3.26) 



and then identify the classical source with the square root of this; quoting again [25], we 
have 

P n(M) = ^^Vdy Y / a t ( 1 )P i (y,7 1 )P i (\x-y\, V ), i = {E, B, EB} , (3.27) 

where the order one Uj coefficients are mildly 7 dependent (within the region I7I < 2 only), 
and we have ignored a bunch of factors 2ir for simplicity. The power spectra appearing in 
this convolution are as usual defined as 



Pb = 
Pe 
Peb 



k 



Pa* 



1 

k 



\A\ 2 , 

A"' 
f 
A 



f 



A*. 



(3.28) 
(3.29) 
(3.30) 



All that matters in this expression is the final time dependence, which can be calculated 
using some rough approximation for the convolution integral for all values of 7 £ [—2,2] (see 
the details in the original papers [25]) and cast in a simple and illuminating form 



Ph(k,Tl) ~ P{k)\T)\ 



2cx 



(3.31) 



with 



This means that 



a(7)=min{4-2| 7 |,3/2}. 



(3.32) 



whence 
S, 



em — 4 



3 

x~ 4 



u' s = - a nu s 



(2 — a + 4ei — 3ae\ — 662)115' 



~j Us ~ y^ cm ~ Pei 



+ 3x 2 Q c 



(3.33) 



(3.34) 



(pp is the energy density of ip). 

3.3 Solving for the Bardeen potentials 

It is now possible to solve the Bardeen equation (3.13); the solutions are obtained as a 
linear combination of two homogeneous solutions, which are ^i(x) = x p J u (x) and ^(^O = 
x p J- l/ (x), with p = 1/2 + C(ei,e2) and v = 1/2 + 0(61,62), and a particular solution 
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of the inhomogeneous equation, obtained through the two homogeneous solutions and the 
Wronskian 

2 



7T 



sin(^7r)a; 2p 1 , 



as 



where 



^ inh (k,x) = «;i(k)*i(s) + w 2 (k)y 2 (x) , 



M = -/ 



dz 



^2(^)^cm(k, Z) 

W(z) 



, 102 (k) = J dz 



fri(z)«S cm (k,z) 
W(z) 



(3.35) 
(3.36) 

(3.37) 



Thence: 



*inh(k, x) = £ dlnz (f ) p [J,(x) ./_„(*) - J„(z) J_„(x)] z 2 5 cm (k, *) 



l r x 

— 2i/ 



1 + 



Av - A Av + 4 



(3.38) 



+ 



ar 



Av + A Av - A 



Z 2 5 em (k, z) , 



plus higher order terms; we keep the second order terms in the expansion since, as we will 
see shortly, the first term does not contribute to the final potentials. Notice here that, unlike 
the convolution Eq. (3.27), where for convenience we rewrote the momentum k in terms of 
the dimensionless x, but the integral was a momentum integral, here we are integrating over 
conformal time rj - this is also why we chose y as dummy variable there, and z here. This time 
integral runs from horizon exit Xh — 1 as this is when the mechanism becomes operational; 
we also have expanded the integrand for large scales x <C 1. 

Focussing only on the "strictly" inhomogeneous solution at z — > x, and ignoring the 
homogeneous solutions as they are subdominant, one can write then 



~ b(k) 



3 — a 



5 — 3a 
3 — a 



Jl5 1 

x 2 p ip (3 - a)(l - a) 



6 - a Pem 
a p^ 



-3(l + 2ei 



n 



s 



x 2 p ip 



(3.39) 
(3.40) 



b{k) denotes the contribution from the scalar field fluctuations at large scales. The last step 
is the quantum-to-classical transition, in which the quantum operators of (3.39) and (3.40) 
are transformed into the squares of their power spectra, leading directly to 



~ b(k) - 



3 — a 
1 



5 — 3a 
3 — a 



(3 - a)(l - a) 



a 



<M*0 



3 (l + 2ei )^fc- 3 / 2 , 



where the dimensionless energy density fraction f2 n is conveniently defined as 

[^n(M)] 2 -^, 

Pip 



(3.41) 
(3.42) 

(3.43) 
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and similarly for 0~. 

Once inflation ends, the Universe enters the reheating phase, which is (hopefully) fol- 
lowed by epochs of matter or radiation domination, once the decay products of the inflaton 
(hopefully) thermalise. The reheating epoch, although now yet completely understood, typi- 
cally lasts a short (on Hubble scales) time, which means that for very large-scale modes which 
are sitting outside the Hubble radius, this process is essentially instantaneous and not influ- 
ential on their dynamics. Hence, in order to tailgate the destiny of the Bardeen potentials to 
our times we need to match the inflationary solutions (3.41) and (3.42) to the corresponding 
results in radiation (or matter) domination 4 . The latter can be found by solving the Bardeen 
equation in such backgrounds 

y» + m > = ^L(U 2 n+)' , (3.44) 

*" + 3n*' = ^(n 2 n+)' , (3.45) 

for radiation and matter domination, respectively, meaning that 

*+(*) = *o + ^ + ^" 3/2 , (3.46) 
*+{x) = * + ^ , (3.47) 



x 3 



in radiation domination, and 



*+(*) = *o + :4 + — a 5 *" 3 ' 2 , (3.48) 



x J ar 



<M*) = * + ^ + ^fc^ 2 , (3-49) 



in matter domination; we have defined = n^fc 3 / 2 / p r>m for radiation and matter, re- 
spectively - note that in radiation domination f]^ oc const whereas in matter domination, 
since Hs oc 1/a 4 and p m oc 1/a 3 then Q^j oc 1/a; useful relations are % = 1/rj in radiation 
domination, and % = 2/rj in matter domination (77 positive). 

Notice the appearance of \I/o and Vl/i, two integration constants which have to be fixed 
by appropriate matching with the preceding inflationary epoch. 

3.4 Matching 

Since the energy-momentum tensor, in the idealised description we adopt, goes through 
a discontinuity when jumping out of inflation, special care must be taken in sewing the 
solutions of the Bardeen equations together. Keeping in mind that we are working in the 
long wavelength limit, the appropriate way to ensure that the gravitational potentials are 
not discontinuous at the sharp transition is to keep the induced three-metric on and extrinsic 
curvature of the spacelike hypersurface which defines the transition continuous. While this 
seems very pompous and technical, in the large-scale limit it simply reduces to demanding 
that the curvature perturbation 

U 2 /*' \ 

)■ (3 50) 



4 We work out the matter domination case for completeness, but this may be a rather academic exercise 
since a very long dust epoch after inflation is hardly ever realised. 
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is conserved through the transition. 

The hypersurface which determines the transition is typically defined by the requirement 
that the energy density on such hypersurface remains constant; we will not dwell into the 
details of this procedure, and simply quote the result 

r=%, (3.5i) 

Pip 

where T defines the transition hypersurface. Eq. (3.51) refers to the transition from inflation 
to matter or radiation, but is applicable to other discontinuities of the same kind. It is 
straightforward to write down the explicit form for T in the matter or radiation cases, which 
looks like 



r 1 



^+2— + 2$ +— ^ 



(3.52) 



where all quantities have to be evaluated at the conformal time corresponding to the transi- 
tion, which we call rfj. 

In order to ensure that the Hubble parameter remains continuous through the transition 
we are automatically led to define the following, self-explanatory, relations: 

V- = ~VT = ~V+ or V- = ~Vt/2 = -»7+/ 2 , (3.53) 

for radiation and matter domination, respectively. 

With the definition of the transition hypersurface, Eq. (3.51) all the necessary ingredi- 
ents to finally pierce together the two halves of the cosmological evolution of the Bardeen 
potentials are in place; in the longitudinal gauge the matching conditions can be cast as 

[T]± = 0, [*]± = 0, [H$ + ^' - {H 2 - H')T}± =0. (3.54) 

Remembering that % 2 — "H' = 2% 2 and "H 2 — "H' = 3"H 2 /2 for radiation and matter domination, 
respectively, one thus obtains the final expressions for the integration constant (we are 
not interested in since it decays rapidly): 

* _ JL { w + 0(a |^ + ^ + + I] n;*-^} ,(3.55) 

(3.56) 

for ultra- relativistic gas, and similarly for dust. The most important contribution comes 
from f2~, while the similar piece is suppressed by powers of the slow roll parameters; 
the largest term Q^j/xj- is not present up to the next order in slow roll, but one can in fact 
show [25] that it vanishes at all orders, as it should. 

In principle one has to ensure that this perturbation does not grow disproportionately, or 
perturbation theory would be invalid. Indeed, while, in this (dis)guise, the constant mode 
seems rather innocuous, the contribution coming from the term could in fact dramatically 
perilous. An easy way to see this is by looking at the relative weight of this perturbed piece, 
which will enter the otherwise null Weyl tensor, with respect to the background, represented 
by the Ricci tensor: 
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it is obvious now that one has to normalise the overall anisotropic contribution in such a 
way that this ratio stays controllably below unity at all times, or at least until today. This 
is automatically true as long as x <C 1 (large scales) and since we also need Q.i{xq-) C 1 or 
the de Sitter background solution does not hold. Notice that, had the 1/xj- mode leaked 
into its contribution would be so much enhanced that, as it is easy to see, it would make 
magnetogenesis and inflation incompatible. 

In any case, this new contribution to the constant mode is to be taken into account when 
inflationary predictions are discussed; in particular, it can lead to non-Gaussian signals, which 
we analyse in our specific model below. 



4 Three-form perturbations 

The outcome of the previous section was that, in what seems to be a fairly general class 
of models of inflationary magnetogenesis, once one proceeds to analyse the first order grav- 
itational perturbations induced by the presence of the magnetic (and electric) fields, there 
is a nonzero transmission coefficient for the magnetic energy density and anisotropic stress 
generated during inflation when moving over to radiation or matter domination. 

We will now repeat the above calculation for the three-form case, pointing out where 
the key differences arise; the non-Gaussian signatures of this scenario will be discussed in the 
next section. 

4.1 Perturbations 

The first step is to write down a Bardeen equation for the three-form system. In order to do 
so, we define the kinetic 

K' = X' + mX] (4.1) 
the background equations for the three-form inflation system are 

x" + mx' + m'x - m 2 x + a 2 v x = o , (4.2) 

* 2 =^G*" 2+ " 2 ^)- < 4 - 3 » 

"'""'"Sil"'^ 1 ' <44) 

In the longitudinal gauge the expressions for the perturbed Einstein tensor of course do not 
change (the three- form is minimally coupled); what does change are the formulas for the 
perturbed energy-momentum tensor. The perturbation equations read 

+ m 2 $ + kH = -ATTGa 2 5p x - 4irGa 2 p cm , (4.5) 

k 2 

V + ZHV' + Wb' + {2W + U 2 )$> -—($-*) = AitGa 2 5P x + AirGa 2 p eUi , (4.6) 

y 

W + %<$> = AirGVxaa - 4vrGa i ^ q emj , (4.7) 
fc 2 ($ - $) = -8^Ga 2 n s . (4.8) 

The three-form perturbations are given by 

B ij = a 3 eij k (i?^ + dkB L ) , Bij k = a^e^ (X + B ) , (4.9) 
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with the transverse (and ultimately only dynamical) degrees of freedom, and where 

k 2 



Vx T 
B' + 3HB + -^aB L 

A 



-aB + K 1 (3*-$) = 0, 



and the energy and pressure density perturbations are given by 



a 2 5p x 



K' 



k 2 

B' + 3HB + -^aB L + K' (3* - $) 



+ a 2 V x (B + 3XV) 



a 2 5P x = -a 2 5p x + a 2 (V XX X + V x ) (B + 3X*) . 



(4.10) 

(4.11) 
(4.12) 



It is straightforward, if not extremely exciting, to combine these equations once more 
into a single second order differential equation for the Bardeen potentials, for instance 



X' 



X 



V,x J 



where the source now is 



3 



A; 2 



S cm = 8nGa 2 [u + [ 2 («' - - (l + ^) 

^Aay VxxX \ ifcJ gem 

We can simplify this again for large scales using 

J( 1 + ^). 2W(ei + 3 , 2) , 

and the sound speed c| = VxxX/Vx to obtain 
d 2 ^ 2(ei + 3e 2 ) d* 



(4.13) 

(4.14) 
.(4.15) 

(4.16) 



dx 2 



x 



dx 



2 2(2ei + 3e 2 ) 
C 5 ^ 



(4.17) 



where the source takes almost the exact same form as in Eq. (3.14); the only relevant difference 
is thus the nontrivial speed of sound, which is in principle time-dependent (although only 
mildly so through the slow roll parameters) : 



2u2 r 



S ~ 



8TrGa z n 



n 



+ 2(l-2 ei -3e 2 )n 5 



A: 2 



-n 5 



k 



(4Pe 



P™) + ^{4 + l)k 2 QM 



18) 



The solution to Eq. (4.17) is again obtainable by resolving the homogeneous equa- 
tion first, and then compose the particular inhomogeneous solution through the Wronskian 
method. The homogeneous solutions are ^i(x) = x p J u (csx) and ^(x) = x p J~v(cs%), with, 
as before, p = 1/2 + 0(e\, e 2 ) and v = 1/2 + 0(e±, e 2 ), and where we have taken the sound 
speed to be constant (up to slow roll corrections); thence, first of all, W(x) = — f sin(i/7r)x 2p_1 
as previously, which then means that the lowest order inhomogeneous solution will not know 
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anything about the nontrivial speed of sound - as is expected since this term is small at large 
scales: 

*inh(k,z) = 2^£ h d]nz(z) p [J v (x)J- v (z)- J„(z)J- v (x)]z 2 S em (Kz) 

=^fidmz(f)*{(f) 



l z 2 S cm (k,z) . 
4.2 Spectra 

While up to this point it looks like little has changed, one has in fact to remember that 
the convolution leading to the source term is very different from the f 2 F 2 case - recall we 
have / = 1 in this example. Indeed, the spectrum of the magnetic field is Pb oc exp(2r), 
see (2.20), for horizon exit IC, and Pb oc exp(2rfc|?7|) for IC at the beginning of inflation, 
Eq. (2.27). The convolution in (3.27) is an integral over momentum; due to the nontrivial k 
dependence of the T coefficient in the exponents, these integrals are very hard to do. However, 
we know that the spectrum is essentially zero except for a narrow band of modes, where we 
can approximate T as a simple power law in k. In the simplest case, when T = const, 
the integral is doable analytically. A difference to the scalar field case is that now there 
is also a contribution from the coupling term to the anisotropic stress in the stress energy 
tensor. However, it is suppressed by the small coupling parameter and we neglect it here for 
simplicity. Furthermore, noticing that Pg = TPeb 3> Peb = FPb 3> Pb we can approximate 
the convolution as 

^ ^ f dyy^-x\ e2T{2 _ y _ lx _ ym ICat _ 1/fc (42Q) 
J a 8 |r/|° 

^ „ f dyy 3 \y - x\ e 2r(y-|y-»|)(i-w»>) , i C at Vi . (4.21) 
J 

The integrals have to be performed only over the range of modes which are amplified, that 
is, up to kA only. Notice also that in the very end what we care about is the final time 
dependence, because this is what is going to tell us how to solve Eq. (4.19). 

The most relevant difference here is that we need to explicitly account for the cutoffs, 
which we call km and kuv-, and the momentum integral will generally depend on these. This 
will change significantly the shape of the result, for no longer the time dependence is going 
to be simplified at the boundaries k = and k = 1/|??|. Therefore, 

Pn ^ M 3 X(fc, k IR , k uv )e 2T ^ 2 - k ^ 2k ^ , IC at - l/k , (4.22) 
P n ^ N 3 X(A;, k IR , k uv )e 2r ^- k ^~^ , IC at m , (4.23) 

where in both cases the functions l(k, kiR, kuv) are polynomials whose shape depends on the 
specific values of the parameters A and k\; their form however does not matter for us since 
they do not depend on time explicitly - moreover, if they are to solve the problem of today's 
observed large-scale magnetisation, they will return enough power at these lengthscales we 
are after. Note that (4.23) is valid as long as T^aM < 1 - were this not satisfied, the 
dominant contribution would be oc I??) 8 . 

We can use these approximated expressions to calculate the form of the source in this 
case. In the language of Eq. (3.31) we would say that a = 3/2 for the two above possibilities, 



1 + 



4z/ - 4 4z/ + 4 



(4.19) 



1 - 



CqZ 



+ 



CqX 



4^ + 4 \v - 4 



and ignoring the mild time dependence coming from the exponential factors - in computing 
Us' from (3.33) this will simply add a subdominant term proportional to Tx to (2 — a). 

We can write down the solutions for the Bardeen potential up to terms of order x°, 
which now contain information about the nontrivial speed of sound: 



~ b(k) - 



3 — a 
1 



1 + 



3a 



3 — a 



n 



X 



(3-a)(l-a) [ n 



+ 2(6 - a)4 - (4 - l)a 3/2 



2a 



(4.24) 



The potentially dangerous mode 0.^/x 2 again cancels away from the final expression for 
\&o just as it should; notice that since there is no sign of cs in the term proportional to 
the anisotropic stress f2^j, then once matched to radiation this contribution is going to be 
suppressed by slow roll parameters. 

One can also look at this solution directly in terms of curvature (, in which case this 
solution reads 



Cinh(£0 ^ 



1 



3 — a ei 



^-3/ 2 _ 2(6-a)4-(4-l)a 3/2 
n 2a p 



(4.25) 



Notice that this contribution is dominant with respect to the "initial" term when the in- 
tegral is evaluated at rji, for in the latter case the exponential factors become relevant and 
rapidly approach zero, killing the contribution (as it must be), whence the goodness of the 
approximated curvature (4.25). 



5 Non-Gaussianities 



The model we have developed introduces a direct coupling between the three-form and the 
one-form, that is, the vector potential. Since in the final expressions for the curvature sourced 
by this coupling (4.25) the vector field always appears quadratically through its energy density 
p cm or the anisotropic stress Us, it is natural to expect some degree of non-Gaussianity in 
the primordial spectrum, even for purely Gaussian vector two-point functions. 

Non-Gaussianity in the context of direct couplings between a scalar or pseudoscalar 
inflaton field and EM has been studied in [20, 27] , from where we will borrow some of the 
machinery for our rough estimations. The relevant difference from their work and ours is the 
form of our coupling term Eq. (2.3), which does not resemble the structure of those most 
popular ones analysed in the literature. There will be some intrinsic non-Gaussianities arising 
from the third order coupled metric and three- form perturbations [46-48] , but we expect this 
contribution to be subdominant, since the direct coupling is strongly enhanced by the rapid 
development of the vector potential instabilities (which in turn give rise to strong magnetic 
fields). 

5.1 The formalism 

Let us now sketch the steps, without repeating all the details which can be easily found in 
the literature, which lead to a reasonable guess of what level of non-Gaussianity we might 
reach in this context. Since the calculations are, euphemistically, tedious, and in view of the 
approximations we are going to employ, we will limit ourselves to one significant example 
wherefrom illustrative conclusions can be drawn. 
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First of all, the leading contribution to the curvature comes from fi~, so what we are 
interested in is the three-point function of this operator, that is (p C m(p)pcm(q)pcm(k)). Each 
p cm at the operator level will contain combinations of |„4(fc)'| 2 , \kA(k)'A(k)\ and k 2 \A(k)\ 2 ; 
we have seen that the parameter space which is relevant for magnetogenesis is that for 
which, schematically: A' ~ TkA > kA since F >1. Thence, we only retain the electric field 
components in T^ u . 

The single a 4 p em (k, i]) will be a convolution of the type 



/ 



d 3 ^ {4(kV A '(|k + k'|)a A (k> A ,(|k + k'\)k'\k + k'\A' (k')A' (\k + k'|) + 3 x h.c.} , 

(5.1) 



A A' 



where by 3 x h.c. we mean the remaining three combinations of creation and annihilation 
operators 5 . The three-point function is going to contain all possible combinations of three 
such terms; let us look at just one of them, that with structure (aaa^aa^a^): 



(PL) Op/ 



d-VdVdV 

(2vr) 9 / 2 



A A' p,p! v v 1 

£i A (k'K A '(|k + k'|) £ f (p'^'dp'-pDerCq'^'dq-q'l) 

(a A (k')o A ,(|k + k , l)a / .(p')<(lp'-p|)4(q')^(|q-q'l)> 

k'\k + k'b'ip' - pk'lq - q'|r(A/)r(|k + k'|)r(p')r(| P ' - P |)r( g ')r(|q - q'|) 

^'^(Ik + k'iK^Mdp'-PlK^Kdq-ql)}, (5-2) 

and a number of other similar terms. We chose this one because it is one nontrivial one since 
it it not simply a product of (p cm ) (/9 2 m ) which would be uninteresting. 

The expectation value (aaa^ aa^ a') returns a series of four combinations of, schematically 
(Spa) 3 x (<5 3 (r — s)) 3 which will lead all to very similar terms; we pick the representative one 



<WV^W<5 3 (P - k'^V - p - q> 3 (k + k' + q - q') 



(5.3) 



This will take care of three summations (over A', //, and v') and two integrals, leaving us with 
one convolution only; the third <5 3 in all terms is nothing else than <5 3 (p + q+k) which enforces 
the "triangularity" of the bispectrum, descending from translation invariance. Finally, we 
deal with the polarisation vectors e with the aid of the formula 

J>?(k)^(k') = l + (^') 2 ' ( 5 - 4 ) 

a/3 



which directly leads to 



1 



So 



d 3 k' 



-k'\k' + k|fc'|k' - p||k' - p||k' - p - q|(5 3 (p + q + k) 



(2vr) 9 / 2 

iWXlk + k'\)r(k')r(\k> - P |)r(|k' - P |)r(|k' - P - q |) 

A(k')A(\k + k'\)A(k')A(\k' - P \)A(\k' - P \)A(\k' - P - q| 
1 + (l + kk'^J 1 + (l - pk'^j 1 + (2 - 2pk' - qk' + pq 



(5.5) 



5 Notice that here we are not taking the expectation value, as it is not a fluctuation, but we retain the full 
vectorial form and bra-ket it only at the end as usual to capture the interferences. 
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where we know that our solutions A(k) are real-valued. 

Some approximations are necessary at this point. We focus on one of the two cases 
examined in Sec. 2, that is, the case where the IC are fixed at rj = —l/k, which means 
that \^2kA — exp T where we neglect its (subdominant for large scales) time dependence. In 
order to facilitate the numerical evaluation of (5.5) we also employ the same simplification we 
adopted for the study of magnetogenesis, and consider T(k) ~ 1/A with a cutoff at k Ak\, 
from which point on T decays, and eventually becomes imaginary at k — > k^. It is a very 
rough approximation, but given the extreme sensibility that we expect from the parameters A 
and k\ in view of the exponential factor, it will serve the scope of this exposition. Collecting 
these choices leads at once to 



(pL) => 



1 



„6r-p6 



d 3 k' 

>:i4a 12 ^ ' J (2vr) 9 / 2 

i+ (i + M') 2 



|k' + k||k' - p||k' - p - q| 5 3 (p + q + k) 



l+(l-j5fc') 2 1 + (^2 - 2pk' - qk' +pq j 



(5.6) 



We study this integral numerically; first of all, we employ the delta function to fix 
k = — (p + q), which eliminates three out of nine variables (excluding the free parameters A 
and k\), and then choose a convenient reference frame - this does not affect the generality 
of the outcome - as 



p=(0,0,p), q= (0,52,93), k= (0,-02, -93- £>)■ 



(5.7) 



At this point since there still are three variables (92, 93, and p), one needs to specify a shape 
for the triangle under inspection. This in general will not guarantee that the resulting non- 
Gaussianity follows this pattern, but one can assess the goodness of the template one tampers 
with. In our case we are essentially interested in looking for possible new shapes, and more 
than anything else we are keen to see whether our model can at once provide magnetogenesis 
and detectable non-Gaussianities. 

One possible way to characterise the level of non-Gaussian signal is of course to look for 
the deviations from simple Gaussian power spectra in the form of its variance. This leads to 
the definition of the /nl parameter for the curvature perturbations £ through its spectrum 
as 



2tt 



(C(p)C(q)) = -^(^(p + q), 
<C(p)C(q)C(k)> = -|(2vr) 5 / 2 /NL^ 2 ( fc )^(P + q + k ) 



10 



p 3 + 9 3 + A: 3 
p 3 q 3 k 3 



£<5 3 (p + q + k). 



Looking at our solution for £ Eq. (4.25) we can infer (£) ~ c(a)(p em ) jt\p^ and obtain 



/nl ^ 



15 



2r-p2i 3 



[c(a)e 2r r 



(Ak A ) 4 



64(2vr) 3 ei 

where we have factorised all shape-related factors in the term 
p 3 + q 3 + A; 3 



^V^(p,q,k) 
PX J 



^(p,q,k) = 



< 3 o 3 /c 3 



I d 3 k'\k' 



+ k||k' - p||k' - p - q| 



(5.8) 

(5.9) 
(5.10) 

L 

(5.11) 
(5.12) 



jl+fl + M') 2 l+[l-pk'^ 2 1 + [2 - 2pk' - qk' + pq^ j ; 
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notice that all momenta appearing in the last Eq. (5.12) are normalised to Ak\, and therefore 
run from to 1. This integral can be fed to numerical solvers once a shape is specified, and is 
expected to not run amok: the all-important amplitude factors are displayed in the prefactor. 

A comment is in order here. In performing the integration we have approximated T as 
a constant in momentum, and cut this at Ak\ because from that point on the actual value 
of T decays as 1/k 2 and is eventually cut off completely shortly thereafter (at k = k\ itself). 
This implies that, since we have six Ts in the integrand, looking at their arguments we see 
that all the momenta involved (the internal k' as well as the external p, q, and k) will need 
to be limited by the same threshold. This implies that, just as it happens for the magnetic 
field which is only generated at the largest scales - for low energy cut off momentum, so 
the non-Gaussianity produced in this model will be limited to large-scale correlations. This 
is a new feature which will help discerning this model from the plenitude of other variants 
available in the literature. 

5.2 Results 

To first determine the shape of the non-Gaussianities at hand, we employ the visualisation 
introduced in [49]. Then we define the quantity 



Without loss of generality, we order the momenta as q < k < p. The relevant function is 
then F(q/p, k/p), which is normalised to unity for the equilateral configuration: F(l, 1) = 1. 
This function quantifies the relative contribution of different non-Gaussian shapes to the 
observable CMB anisotropies. To avoid plotting the same region twice, we plot only the 
region q < k, q > 1 — k. The result is shown in Figure 5.2. 

We see that the non-Gaussianity peaks at the equilateral shapes. This is what we 
might have expected, since equilaterality is the typical outcome for perturbations generated 
within the horizon: in the present model we have seen that since the exponential terms, 
responsible for the perturbation generation, blow up quickly at the beginning after we set the 
IC and then become constant. Having verified this numerically, we will restrict to equilateral 
configurations in the following and set q = k = p. 




(5.13) 



Scale dependence 



F(k,q) 




k 



1.0 




0.0 



10 4 



100 



10' 



0.001 



0.01 



0.1 



10 



100 



q 



o.o 



p 



Figure 3. Right. The non-Gaussian shape function F(q, k) defined in the text. In our case it 
clearly peaks at the equilateral configuration q = k = p. Left. The scale dependence for equilateral 
non-Gaussianity. 



In order to grasp what the scale dependence is going to be, we plot the an arbitrarily 
normalised non-Gaussianity in the left panel of figure 5.2. For large scales, Jnl remains a 
small constant, but for the subhorizon scales there is a strong amplification with the increasing 
wavemode. The numerical value of non-Gaussianity will depend upon the model parameters, 
and in the rest of this section we study the compatibility of the model parameters with the 
two-fold requirement of both sufficient magnetic field generation and non-Gaussianity that 
lies in the interesting range of roughly 1 < |/jvx| < 100. 

In the specific case under scrutiny said parameters are the UV cutoff scale k\ and the 
strength of the coupling A. As is expected (recall the study of magnetogenesis in Sec. 2 the 
result will be exponentially sensitive to variations in A; in waiting for the results of the full 
numerical integration (also with the numerical solution of the Euler-Lagrange equation (2.15) 
- recall we are working with an idealised de Sitter space) what we want to show here is that 
it is possible to generate large-scale magnetic fields and, at the same time, realise a detectable 
beyond-Gaussian imprint in the CMB. 



log /nl prefactor log 8b ® 1 / Mpc 




-2.10 - 2.05 - 2.00 -1.95 -1.90 -2.04 -2.02 -2.00 -1.98 -1.96 



log A log A 

Figure 4. Left. We zoom onto the most interesting region of parameter space of our Fig. 1 (same 
colours), while we superimpose the values of the prefactor to /nl from Eq. (5.12) in a logarithmic 
scale. Right. Here make contour plots for the regions where this coefficient falls in the ranges: smaller 
than 10~ 2 (light blue, uninteresting), 10~ 2 to 1 (yellow, small), 1 to 100 (green, detectable), and 
beyond 100 (light red, dangerous); on top of these regions are drawn the curves of constant magnetic 
field strength at 1/Mpc (in Gauss units). 

In Fig 4, left, we zoom onto the most interesting region of the (A, k\) parameter space, 
that where the generated magnetic fields could match those required by observations on 
large scales. The light blue region is that for which the field strength is below 10 -25 Gauss at 
1/Mpc, and is uninteresting phenomenologically; the yellow region goes up to 10" 15 Gauss, 
and would, with the help of some dynamo mechanism, make up for observations; the green 
region has a field stronger than 10 -15 and can directly account for Mpc scale magnetisation; 
the light red region is excluded due to background backreactions constraints. The coefficient 
f° r /nl which can be read off Eq (5.12) can vary significantly in this region, due to the 
exponential sensitivity to the coupling and the high power with which the cutoff scale appears; 
this coefficient then will be enhanced (suppressed) by up (down) to a couple of orders of 
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magnitude when the full integration and momentum dependence is considered, but it still is 
useful to understand the basic picture. We see how it is possible to obtain order one non- 
Gaussianities in the green and yellow areas: this was highly nontrivial in the beginning; we 
plan to explore the details of such signatures (shapes, different potentials, realistic inflationary 
set-up, more detailed backreaction bounds, tensor perturbations) in a forthcoming work. 

The right hand side of Fig 4 instead shows the regions where the /nl coefficient is 
below 0.01 (light blue, deemed uninteresting), between 0.01 and up to 1 (yellow, small non- 
Gaussianity) , between 1 and 100 (green, detectable non-Gaussianity), and beyond 100 (light 
red, too much power in the bispectrum). The yellow and green bands are going to be shifted 
to the left or right, but only by some very small amount since we do not expect the integral J* 
to be much different from one. On top of these regions we show the values for the generated 
magnetic field at the Mpc scale. This is in some sense the inverse or complementary plot 
to the left hand side one, and once more demonstrates how we can find some bands where 
the non-Gaussian signal is strong enough to be measurable, and at the same time large scale 
magnetisation is successfully achieved. 

What we read off these figures is the effect of the natural cutoff of three-form inflation. 
In the first place this affords to amplify extremely rapidly the EM field once the univocal 
U (l)-invariant coupling is introduced, without incurring in the most basic backreaction effect 
of blocking de Sitter expansion. Since the scale k\ can be chosen to lie in the low-energy 
range, the band of magnetic field which is boosted is very narrow; the spectrum also decays 
in the IR, ensuring that no divergence is hiding there for any duration of inflation. 

Secondly, the UV threshold limits from above the size of the non-Gaussian scale which 
is going to be processed by the coupling. This limit automatically gives birth to interesting 
phenomenology, shapes that are not obvious, and peculiar running with the scale, mimicking 
that of the generated EM field. It is instructive to compare our findings with those of [20]: 
in their case the logarithmic enhancement of the bispectrum is due to a very long period of 
inflation before the CMB scale has exited the horizon, and is a result of the vector waves 
being amplified from very low energies (largest scales) which all contribute to the bispectrum 
since their spectrum is flat. In our scenario on the other hand, for k\ not tremendously larger 
than the Hubble scale today, everything happens only at the scales of interest. In the end it 
seems possible to reconcile magnetogenesis and its epiphenomenal non-Gaussianity. 

6 Summary 

Three-form inflation has proven to be an enticing arena for the phenomenology of EM cou- 
plings to inflation. We have analysed in more detail the unique gauge-invariant coupling 
which is allowed between the three-form and the EM one-form, going beyond our previous 
background findings. The outcomes are bold and encouraging us to undertake further work 
in this direction. 

First of all, "three-magnetogenesis" proves to be perhaps the only one magnetogenesis 
model which eschews strong coupling and background backreaction bounds, yet is capable 
of amplifying EM fluctuations on large scales to a degree which would be compatible with 
empirical evidence. 

Secondly, at the first-order perturbation level one finds, in agreement with recent works, 
that the constant mode of the Bardccn potential and the curvature perturbations after infla- 
tion ends are influenced directly by the spectrum of the produced vector fluctuations on large 
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scales. In the simplest case we have described such fluctuations are statistically isotropic, 
but one can easily imagine extensions to more general cases. 

In the third place, the "revised" curvature perturbation, containing bilinears in the 
Gaussian EM spectra, is genuinely non-Gaussian. We have explored in some detail the 
simplest model to find captivating hints to nontrivial new signatures, shapes, and amplitudes 
within the detectable span. What we found was scale-dependent equilateral non-Gaussianity 
that is within the viable and future-observable range, where at the same time the desired 
amount of magnetic field is generated. Even though more precise predictions are demanded 
to be of use for data analysis, the toy-model we have adopted for the latter part of this 
study already contains all the important physics, in particular the central role played by 
the three-form duality, which translates in a low-energy UV cutoff of the enhanced vector 
spectra. 

This model's features permit a direct link between magnetogenesis, inflation, and bis- 
pectra, in a way that can be reconciled with multiple observational pieces of evidence, a task 
definitely not paltry; the positive results of our analysis prompt us to go beyond the na ive an- 
alytical estimations we have presented, in order to fully eviscerate the precise phenomenology 
of this model. 
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